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Abstract 

The reaction e + e~ — > e + e~pp is studied with the L3 detector at LEP. The 
analysis is based on data collected at e + e~ center-of-mass energies from 183 GeV 
to 209 GeV, corresponding to an integrated luminosity of 667 pb~ x . The 77 — > pp 
differential cross section is measured in the range of the two-photon center-of-mass 
energy from 2.1 GeV to 4.5 GeV. The results are compared to the predictions of 
the three-quark and quark-diquark models. 
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1 Introduction 



Electron-positron colliders are a suitable place for the study of two-photon interactions via the 
process e + e~ — > e + e~7*7* — > e + e~X, where 7* denotes a virtual photon. The outgoing electron 
and positron carry almost the full beam energy and are usually undetected, due to their small 
transverse momenta. The final state X has, therefore, a low mass as compared to the e + e~ 
center-of-mass energy, y/s, and its transverse momentum is almost zero. The small photon 
virtuality allows the extraction of the cross section 17(77 — y X) for real photon collisions, once 
the photon flux is calculated by QED [Tj. 

Calculations of the cross section 77 — > baryon antibaryon were performed using the hard 
scattering approach of Brodsky and Lepage [2]. In this formalism, the process is factorized 
into a perturbative 77 — > qq amplitude and a non-perturbative part described by the quark 
distribution functions of the baryon. Such calculations with three-quark distribution functions 
[H1H| yielded results about one order of magnitude below early 77 — > pp measurements [5] for 
77 center-of-mass energies, W^, less than 3 GeV. The quark-diquark model was proposed 
as a possible way to model non-perturbative effects. Within this model, the partonic structure 
of the baryon is described by a quark-diquark system rather than three quarks. The composite 
nature of diquarks is taken into account by form factors ensuring that the three-quark model 
is recovered at an asymptotically large momentum transfer. 

This Letter presents a study of the differential cross section of the reaction 77 — ► pp in the 
W 77 range from 2.1 GeV to 4.5 GeV. The data sample corresponds to a total integrated lumi- 
nosity of 667 pb _1 collected with the L3 detector [Zj at y/s = 183 — 209 GeV. The analysis is 
based on the central tracking system and the high resolution BGO electromagnetic calorimeter. 
The events are selected by the track triggers |B]. 

The pp final state in two-photon collisions was previously studied at lower y/s 0E1 and 
recently at LEP [TU] with lower statistics. The present study extends the W 77 range. Our results 
are compared to these experiments and to recent theoretical predictions of the quark-diquark 
model [TT] . 

2 Monte Carlo simulation 

Monte Carlo events for the 77 — > pp reaction, as well as for the background processes 77 — > K + K 
and 77 — ► 7r + 7r~, are generated with the PC Monte Carlo program ^2], for each beam en- 
ergy, within the formalism of Reference 03 The e + e~ — > e + e~r + r _ , e + e~ — > e + e~/i + /z~ and 
e + e~ — > e + e~e + e~ background processes are simulated with the DIAG36 ^2] generator which 
includes the full set of 0(a 4 ) QED diagrams. The generated events are passed through the 
full L3 detector simulation using the GEANT ^3] and GEISHA j^] programs and are re- 
constructed with the same programs as the data. Time dependent detector inefficiencies, as 
monitored during the data taking period, are taken into account. 

3 Proton-antiproton event selection 

In order to match the experimental acceptance with the range of the theoretical predictions, the 
data are analysed only in the region | cos#*| < 0.6, where 8* is the center-of-mass production 
angle of the proton. Events are first selected by requiring two well-reconstructed tracks of 
opposite charge. The track selection criteria are: 
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a distance of closest approach to the interaction point less than 3 mm in the plane 
transverse to the beam direction, 



• at least 30 hits, out of a maximum of 62, in the tracking chamber, 

• a matched energy cluster in the electromagnetic calorimeter, 

• a transverse momentum, p t , greater than 400 MeV, to ensure a high trigger efficiency and 
electron rejection. 

The identification of 77 — > pp events is mainly based on three artificial neural networks, 
used to separate antiprotons from e~, \x~ and h~, where h~ represents either a 7r~ or K~. Each 
neural network consists of five input nodes, a single layer of five hidden neurons and one output 
neuron. The following measured quantities are associated with the five input nodes: 

• the momentum of the antiproton, 

• the probability for the proton mass hypothesis based on the mean energy loss dE/dx 
measured in the tracker, 

• the ratio E t /pt, where E t is the transverse energy measured in the electromagnetic 
calorimeter, 

• the number of BGO crystals in the calorimetric cluster associated with the antiproton, 

• the ratio between the energy deposited in the central crystal of the cluster and the sum 
of the energies deposited in the three- by-three matrix of crystals around it. 

The two last variables exploit the typical signature of an antiproton annihilating in the BGO 
electromagnetic calorimeter: a broad shower spanning several crystals, as opposed to narrow 
showers for electrons, minimum ionisation deposition for muons and compact showers for low 
energy hadrons. Each neural network is trained with the corresponding sample from Monte 
Carlo simulations so that its output value is close to one for antiprotons. 

A particle is identified as an antiproton after a cut on the output of the three neural 
networks. The efficiency for antiproton detection is found to be 74%. Even though the e~ 
neural network rejects more than 98% of the electrons, the e + e~ — > e + e~e + e~ cross section 
is more than two orders of magnitude greater than the e + e~ — > e + e~pp cross section and the 
electron contamination still remains important. To reduce this background, the ratio E t /p t of 
the proton candidate is required to be less than 0.6, as shown in Figure ^ In addition, the 
confidence level that the dE/dx measurement is consistent with a proton must be more than 
5%. These cuts eliminate more than 95% of the remaining e + e~ — > e + e~e + e~ background. 

4 Exclusive 77 — > pp event selection 

To select exclusive 77 — > pp events the following further cuts are applied: 

• Events with a photon candidate outside a cone of 36 degrees half-opening angle with 
respect to the antiproton direction are rejected. A photon candidate is defined as a 
shower in the electromagnetic calorimeter with at least two adjacent crystals, an energy 
greater than 100 MeV and no charged tracks within a cone of 200 mrad half-opening 
angle. 
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• The square of the total transverse momentum of the proton-antiproton pair (YPt) 2 
must be less than 0.1 GeV 2 . This cut reduces inclusive background events of the type 
e + e~ — > e + e~ppX, where X represents one or more unobserved particles. 

These cuts yield a total number of 989 selected events. Figure 121 shows the distribution of 
the effective mass of the pp pair, identified with W 77 . It is obtained by assigning the proton 
mass to the two tracks. 

The exclusive background is estimated by processing the corresponding Monte Carlo events 
with the same analysis cuts. The 77 — > 7r + 7r~ and 77 — > K + K~ cross sections given in Reference 
IT61 are used. This background is found to be negligible in the region W 77 < 2.5 GeV and 
increases up to 25% above 2.5 GeV. The e + e~ — > e + e~e + e~ background varies from 5% in the 
low mass region to 30% in the high mass region. The e + e~ — > e + e~/i + /i~ contamination is less 
than 0.2% and neglected. The contamination of inclusive channels and of e + e~ — > e + e~r + T~ 
is estimated by fitting the tail of the (YPt) 2 distribution. The extrapolation to the region 
(J2pt) 2 < 0.1 GeV 2 gives an average background level of (2 ± 1)%. The sum of all backgrounds 
is subtracted bin by bin for the cross section determination. The background composition of 
each bin is detailed in Tabled 

5 Cross section measurements 

The detection efficiency is determined by Monte Carlo simulation in bins of W 77 and | cos 9* \ . It 
takes into account the track acceptance and selection criteria, the exclusive proton-antiproton 
identification criteria and the track trigger efficiency This trigger efficiency is 91% at W 77 = 2.1 
GeV and rises to 98% for W 77 > 3 GeV. The higher level trigger efficiencies are estimated from 
the data themselves, using prescaled events, and range from 88% at W 77 = 2.1 GeV up to 
99% for W^j > 3 GeV. The total selection efficiency is found to be maximum, about 6%, at 
~ 2.5 GeV and | cos^*! < 0.1. It decreases both at larger | cos6**| and W 77 , due respectively 
to limited angular acceptance and antiproton identification efficiency. The average efficiencies 
as a function of for | cos#*| < 0.6 and as a function of | cos#*| for PV 77 < 4.5 GeV are 
presented in Figure El 

The cross sections are evaluated in bins of W 77 and | cos 9* | for different values of yfs. Due 
to the limited selection efficiency near threshold (Figure EJ, the measurement is restricted to 
the range W 77 > 2.1 GeV. A total of 938 events are selected after this cut. The differential cross 
sections are integrated to obtain the production cross sections e + e~ — > e + e~pp for | cos#*| < 0.6 
and 2.1 GeV < W 77 < 4.5 GeV as a function of y/s. Since the results show no significant y/s 
dependence they are combined into a single measurement at (a/5) = 197 GeV: 

a(e + e" -> e + e"pp) = 26.7 ± 0.9 ± 2.7 pb 

where the first uncertainty is statistical, the second systematic. 

The following systematic uncertainties are considered. The uncertainty due to the selection 
procedure is evaluated by varying the selection cuts. The uncertainty on the neural network 
outputs is estimated by varying the input variables according to their resolution. An uncertainty 
of 50% on the 77 — ► 7r + 7r~ and 77 — > K + K~ cross sections is propagated in the background 
subtraction. The uncertainty on the e + e~ — > e + e~e + e~ background is estimated to be 30%. 
The uncertainty of 50% on the level of inclusive background is also included. For low values of 
iy 77 the dominant systematic uncertainty is about 7% and is due to selection cut variation. It 



4 



remains the most important systematic uncertainty for W yi < 3.2 GeV, as large as 30%. For 
Wj-y > 3.2 GeV, the background subtraction uncertainty becomes the most important, rising 
to 55% in the last JV 77 bin. The systematic uncertainties shown in Table El are the quadratic 
sum of the different sources. 

The differential cross section da( , yy — > pp)/d|cos#*| for real photon-photon collisions is 
extracted as a function of W 77 by dividing out the two-photon luminosity function and extrap- 
olating to Q 2 = with a GVDM form factor . The luminosity functions are evaluated for 
each -y/i and W 77 interval. The measured differential cross sections are integrated to obtain the 
cross sections a( r y r y — > pp). An additional uncertainty of 5%, due to the choice of the photon 
form factor, is included in the systematics. The results are reported in Table 121 The present 
measurement is of higher statistical precision and extends towards higher values of than 
the previous results [3101110] • Agreement is observed within the quoted uncertainties, except 
for the measurements of Reference EH which lie below our data in the low mass region. 

6 Discussion of the results 

The results are compared to the predictions of the three-quark [Hj and the recent quark-diquark 
model |Tj] in Figure EJ The three-quark prediction is based on the leading order QCD calcu- 
lations, using the distribution function of Chernyak and Zhitnisky JH]- The quark-diquark 
calculation is performed with the standard proton distribution amplitudes JT] and includes 
first order corrections due to the non-vanishing proton mass. While the shapes of the theo- 
retical curves are quite similar, the normalisations are significantly different. The predictions 
of the three-quark model are about an order of magnitude below the measurement, whereas 
the quark-diquark predictions describe the data much better. The apparent change in the 
logarithmic slope of the cross section observed in data at W 77 around 3 GeV is, however, not 
reproduced by this model. 

In order to investigate this further, the differential cross sections are summed separately 
in three mass intervals: 2.1 GeV < W yy < 2.5 GeV, 2.5 GeV < W 77 < 3.0 GeV and 3.0 GeV 
< iy 77 < 4.5 GeV. The results are reported in Table HU and plotted in Figure El A distinctive 
difference between the three distributions can be observed. No prediction is available for the 
diquark model for W 77 < 2.5 GeV, but it can be seen that the data has a qualitatively different 
behaviour to the diquark predictions, as it is strongly peaked at large angles. Assuming the 
presence of a single angular momentum state in the s-channel of 77 — > pp production, a 
fit to a single spherical harmonic: V °, , ^2 or ^2 * s performed. This does not give a 
satisfactory result. An acceptable fit is obtained using all of them, but only the Y 2 ° contribution 
is significantly different from zero. A satisfactory fit is also obtained using only the V ° and V 2 ° 
harmonics with the fractions 8% and 92% respectively, as shown in FigureEK- The intermediate 
region exhibits a rather flat dependence, which partially agrees with the model predictions. 
The forward peaking behaviour of the differential cross section in the high mass interval is well 
reproduced by the quark-diquark model. Only in this region, then, the data can be described 
by the Brodsky-Lepage hard scattering approach. 

The presence of two distinct production mechanisms can also be seen in Figure and 
Table 0] by considering separately the cross sections as a function of in a large angle 
region, | cos#*| < 0.3, and in a region 0.3 < | cos#*| < 0.6. The shape of the efficiency as a 
function of is similar for the two regions. The large angle cross section dominates in the 
low mass region and shows a steeper fall-off with W / 77 than the small angle cross section. A fit 
of the form a(W) oc W~ n to the small angle region shown in Figure |Hb gives a reasonably good 
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description with n = 9.8 ± 0.3. The data are also in agreement with the quark-diquark model 
prediction in this region. On the other hand, Figure Efc, the large angle cross section cannot be 
described by a simple power law behaviour for W 77 less than 3 GeV and does not follow the 
quark-diquark prediction. In this case, the change of shape at 3 GeV is similar to, but more 
pronounced than that in Figure 0] 

In conclusion, an accurate study of the 77 — > pp process is performed. The present results 
are of higher statistical precision and extend towards higher values of W 77 than those of previous 
experiments. Current models are only moderately successful in describing the observed data. 
These data can provide useful inputs for predictions of other 77 — > baryon antibaryon channels 

EH- 
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Table 1: The estimated inclusive, e + e and h + h backgrounds in each W 77 bin. 



w 77 


(Wrr) 


Number of 


Background 


<7( 77 _> pp) 


(GeV) 


(GeV) 


Events 


(%) 


(nb) 


2.1 - 2.2 


2.15 


216 


4 


5.35 ± 0.36 ± 0.55 


2.2 - 2.3 


2.25 


252 


5 


4.34 ± 0.27 ± 0.41 


2.3 - 2.4 


2.35 


182 


7 


2.86 ± 0.21 ± 0.28 


2.4 - 2.5 


2.45 


111 


10 


1.78 ± 0.17 ± 0.18 


2.5 - 2.6 


2.55 


61 


16 


1.01 ± 0.13 ± 0.11 


2.6 - 2.8 


2.69 


57 


25 


0.50 ± 0.07 ± 0.09 


2.8 - 3.2 


2.97 


32 


54 


0.12 ± 0.02 ± 0.05 


3.2 - 3.6 


3.37 


15 


57 


0.06 ± 0.02 ± 0.03 


3.6 - 4.5 


3.95 


12 


69 


0.02 ± 0.01 ± 0.01 



Table 2: The number of events, estimated background and 77 — > pp cross section as a function 
of Wjj for I cos 9* I < 0.6. The average value (W 77 ) of each bin corresponds to a weighted average 
according to a W~J- 2 distribution. The first uncertainty is statistical, the second systematic. 



cos#*| 



da(77 — > pp)/d| cos#*| (nb) 



2.1 GeV < W~ n < 2.5 GeV 



2.5 GeV < W 77 < 3.0 GeV 



3.0 GeV < W~ n < 4.5 GeV 






- 0.1 


0.1 


- 0.2 


0.2 


- 0.3 


0.3 


- 0.4 


0.4 


- 0.5 


0.5 


- 0.6 



10.36 ± 0.64 ± 1.08 
7.79 ± 0.57 ± 0.81 
6.15 ± 0.53 ± 0.64 
5.39 ± 0.54 ± 0.56 
3.77 ± 0.50 ± 0.39 
2.30 ± 0.48 ± 0.24 



0.86 ± 0.16 ± 0.14 
0.56 ± 0.12 ± 0.09 
0.62 ± 0.14 ± 0.10 
1.01 ± 0.18 ± 0.17 
0.82 ± 0.18 ± 0.14 
1.03 ± 0.22 ± 0.17 



0.06 ± 0.02 ± 0.03 
0.04 ± 0.02 ± 0.02 
0.03 ± 0.02 ± 0.02 
0.06 ± 0.02 ± 0.03 
0.07 ± 0.03 ± 0.04 
0.11 ± 0.04 ± 0.05 



Table 3: The differential cross section as a function of | cos#*| for the different W iy ranges. The 
first uncertainty is statistical, the second systematic. 
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W 77 (GeV) 



^(77 -» PP) (nb) 



cos#* < 0.3 



0.3 < cos0* < 0.6 



2.1 
2.2 
2.3 
2.4 
2.5 
2.6 
2.8 
3.2 
3.6 



2.2 
2.3 
2.4 
2.5 
2.6 
2.8 
3.2 
3.6 
4.5 



2.15 
2.25 
2.35 
2.45 
2.55 
2.69 
2.97 
3.37 
3.95 



3.58 


± 


0.27 ± 


2.95 


± 


0.21 ± 


2.01 


1 

± 


0.17 ± 


1.21 


1 

± 


0.13 ± 


0.54 


1 

± 


0.09 ± 


0.19 


± 


0.04 ± 


0.03 


± 


0.01 ± 


0.02 


± 


0.01 ± 


10 


± 


4 ± 



6) xlO" 



1.78 ± 0.27 ± 0.18) 
1.40 ± 0.18 ± 0.13) 
0.85 ± 0.13 ± 0.08) 
0.57 ± 0.11 ± 0.06) 
0.47 ± 0.09 ± 0.05) 
0.31 ± 0.06 ± 0.05) 
0.09 ± 0.02 ± 0.04) 
0.04 ± 0.01 ± 0.02) 
11 ± 5 ± 7) xlO 



-3 



Table 4: 



77 



pp cross section as a function of W 77 for | cos^*| < 0.3 and 0.3 < | cos 9*\ < 0.6. 



The average value (W 77 ) of each bin corresponds to a weighted average according to a W, 
distribution. The first uncertainty is statistical, the second systematic. 
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Figure 1: The ratio between the transverse energy deposited in the electromagnetic calorimeter, 
E t , and the transverse momentum, p t , for the proton candidate after the antiproton selection. 
The e + e~ — > e + e~e + e~ Monte Carlo prediction is superimposed on the data. 




W w (GeV) 



Figure 2: The effective mass of the pp pair, W 77 , for the 989 selected events. The 938 events 
at the right of the cut are used in the subsequent analysis. 
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Figure 3: The selection efficiency as a function of a) for | cos#*| < 0.6 and b) cos#* (full 
circles) for 2.1 GeV < W 77 < 4.5 GeV. The full squares indicate the efficiency for detecting 
two charged tracks of opposite charge in the detector, with no further requirements. 
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Figure 4: The 77 — ► pp cross sections as a function of W 71 for | cos6>*| < 0.6 compared to 
the three-quark model calculation e 3 and to the recent quark-diquark model prediction [11J, 
available for W 77 > 2.5 GeV. Statistical and systematic uncertainties are added in quadrature. 
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